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Abstract 

The  transport  properties  of  the  membrane  and  the  related  dissipative  phenomena  significantly  affect  the  performance  of  polymeric  electrolyte 
membrane  fuel  cells.  In  this  work  the  transport  mechanisms  taking  place  inside  the  membrane  and  the  interaction  of  these  with  the  adjacent 
electrodes  will  be  analysed  in  some  detail  to  gain  a  better  understanding  and  control  of  cell  performance  and  particularly  of  open-circuit  voltage 
losses,  which  can  be  anything  up  to  30%  of  the  effective  voltage  of  an  operating  cell,  as  such  an  understanding  could  lead  to  significant  technological 
improvements. 

Here  we  present  an  analytical  and  interpretative  description  to  try  to  identify  the  irreversible  electrode  phenomena  under  open-circuit  conditions 
in  terms  of  a  dissipative  steady-state  and  related  “parasitic  currents”. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

It  is  well  known  that  ionic  conduction  in  polymeric  proton 
conductive  materials  such  as  Nation®  takes  place  under  condi¬ 
tions  which  are  significantly  different  from  those  occurring  in 
an  electrolytic  solution  [1-3]. 

•  First  of  all,  only  the  positive  charges  are  able  to  move  and  to 
carry  current,  while  the  negative  ones,  being  an  integral  part 
of  the  solid  structure  of  the  membrane,  cannot  move  at  all. 

•  Secondly,  the  water-drag  phenomena  associated  with  positive 
charge  migration  are  more  important,  the  flow  ratio  between 
protons  and  water  molecules  being  about  2.3.  Positive  charges 
do  not  migrate  as  free  protons,  but  as  solvated  species  and  each 
positive  ion  consists  of  one  proton  and  two  or  three  molecules 
of  solvent. 

•  Thirdly,  as  a  consequence  of  the  water-drag  phenomena,  the 
overall  water  flux  through  the  membrane  consists  of  two  or 
more  contributions:  one  due  to  solvation  (“drag”  term),  which 
is  directly  associated  with  proton  migration  and  the  electric 


*  Corresponding  author.  Tel.:  +39  010  3532926;  fax:  +39  010  3532589. 
E-mail  address:  betta@diam.unige.it  (E.  Arato). 

0378-7753/$  -  see  front  matter  ©  2006  Elsevier  B.V.  All  rights  reserved, 
doi :  10. 1016/j  .jpowsour. 2005 . 1 2.093 


current,  and  another  or  others  connected  with  non-uniform 
membrane  hydration  (“diffusive”  term)  or  to  the  pressure  gra¬ 
dient  (“forced”  term). 

•  Fourthly,  water  cannot  simply  be  considered  as  an  inert  sol¬ 
vent;  on  the  contrary,  the  hydration  level  of  the  membrane  and 
the  water  fluxes  affecting  hydration  are  important  variables  on 
which  the  transport  characteristics  of  the  membrane  such  as 
conductivity,  diffusivity  and  permeability  strongly  depend. 

•  Finally,  an  apparently  different,  but  probably  connected,  and 
perhaps  even  more  intriguing  question  regards  the  open- 
circuit  voltage  losses  (OCVL)  of  a  polymeric  electrolyte 
membrane  fuel  cell  (PEMFC).  The  observed  open-circuit 
voltages  are  about  0.2  V  less  than  the  ones  expected  according 
to  equilibrium  conditions  (Nernst  equation).  These  OCVL  are 
generally  and  vaguely  attributed  to  “dissipative  phenomena”, 
which,  obviously,  also  have  to  act,  and  in  particular,  under 
open-circuit  conditions.  On  the  contrary,  when  an  interpreta¬ 
tion  in  terms  of  “parasitic  currents”  is  invoked  [9,10],  little 
effort  is  really  made  to  explain  how  these  currents  can  occur 
and  act  in  an  open  circuit. 

The  further  decrease  in  the  cell  voltage  under  closed-circuit 
conditions  is  better  understood:  the  related  dissipative  phenom¬ 
ena  have  been  much  more  closely  studied  and  clearly  determined 
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Nomenclature 

b 

see  (A. 4) 

C 

concentration,  referring  to  the  aqueous  phase 

(kmolm  3) 

C' 

concentration,  referring  to  the  membrane 

(kmol  m-3) 

cl 

a-dimensional  diffusivity,  see  (9) 

D 

diffusivity  (m2  s-1) 

e 

swelling  parameter  of  the  membrane 

F 

Faraday’s  constant  (A  s  kmol- 1 ) 

I 

specific  current  (A  m-2) 

k 

equilibrium  coefficient,  see  (A.l) 

n 

a-dimensional  flux,  see  (9) 

N 

number  of  moles  (kmol) 

N 

flux,  referring  to  the  aqueous 
(kmol  m-2  s-1) 

phase 

N' 

flux,  referring  to  the  membrane  (kmol  m- 

2s-1) 

P 

absolute  pressure  (J  m-3) 

R 

gas  constant  (J  kmol- 1  K- 1 ) 

T 

absolute  temperature  (K) 

V 

voltage  (J  A-1) 

V 

volume  (m3) 

z 

spatial  co-ordinate  (m) 

Greek  letters 

a 

mean  degree  of  solvation 

p 

transfer  coefficient 

£ 

volumetric  fraction  of  the  aqueous  phase 

a-dimensional  voltage,  see  (9) 

<p 

volumetric  fraction,  see  (2) 

X 

degree  of  hydration,  see  (2) 

71 

a-dimensional  pressure,  see  (9) 

$ 

a-dimensional  co-ordinate,  see  (9) 

Subscripts 

A 

anion 

c 

diffusive 

i 

generic  component 

0 

open  circuit 

r 

reference 

t 

total  water 

w 

water 

2 

bisolvated  proton 

3 

trisolvated  proton 

a 

solvated  proton  mean 

ti 

electrical 

X 

hydration 

71 

pressure 

in  terms  of  concentration  polarisation,  activation  polarisation 
and  electrical  resistances  of  the  membrane.  Further  considera¬ 
tions  on  the  mass  transfer  inside  the  gaseous  phases,  which  is 
responsible  for  the  concentration  polarisations  and  current  lim¬ 
itations,  will  be  specifically  discussed  elsewhere  [4,5].  Here  the 


problem  of  the  membrane  resistance  will  be  discussed  in  some 
detail,  in  particular  from  a  rather  fundamental  point  of  view. 
For  instance,  the  explicit  consideration  of  the  migration  of  the 
solvated  molecules  (H20)2H+  and  (H20)3H+  could  throw  more 
light  on  the  transport  properties  of  the  membrane  and  the  related 
dissipative  phenomena.  In  brief,  it  is  hoped  that  any  clarification 
of  the  transport  mechanisms  taking  place  inside  the  membrane 
and  their  interactions  with  the  adjacent  electrodes  will  be  use¬ 
ful  for  understanding  and,  then,  controlling  the  performance  of 
PEMFCs  as  they  become  increasingly  commercialised.  In  par¬ 
ticular,  any  path  towards  understanding  OCVL,  which  can  be 
anything  up  to  30%  of  the  effective  voltage  of  an  operating 
cell,  can  be  an  important  premise  for  significant  technological 
improvements. 

2.  Mass  and  charge  transfer  through  the  membrane 

The  membrane  is  considered  here  as  a  one-dimensional  het¬ 
erogeneous  system  made  of  a  solid  and  an  aqueous  phase.  The 
z  co-ordinate  is  oriented  through  the  membrane  thickness  from 
the  anode  (2  =  0)  to  the  cathode  (z  =  s).  In  the  aqueous  phase,  the 
following  components  are  considered: 


•  water  H2O  (index  w); 

•  the  fixed  negative  ions  (index  A); 

•  the  bisolvated  proton  (H20)2H+  (index  2); 

•  the  trisolvated  proton  (H20)3H+  (index  3); 

•  as  an  alternative  to  the  last  two  components,  the  mean  solvated 
proton  (H20)„H+  (index  a). 


On  the  contrary,  free  H+  protons,  mono-solvated  (H20)H+ 
protons  and  multi-solvated  protons  are  considered  less  important 
and  are  neglected.  The  dissolved  neutral  components  such  as  H2 
and  O2  are  similarly  neglected  initially  because  of  their  very  low 
solubility. 

The  concentrations  Ca  (or  C2,  C3),  Cw,  Ca  refer  to  the  unit 
volume  of  the  aqueous  phase  inside  the  pores  of  the  mem¬ 
brane,  while  the  total  equivalent  concentration  of  this  aqueous 
phase,  Cw,  can  be  considered  approximately  constant  and  cor¬ 
responding  to  the  inverse  of  the  molar  volume  of  liquid  water 
(Cw  ^  1/18  mol  cm-3).  In  such  a  way  swelling  can  be  assumed 
to  be  proportional  to  the  total  hydration  of  the  membrane. 


AV  = 


fl  wt 

Cwt 


Mi  _  A/;Cwt 
AV  «wt 


(1) 


In  a  simplified  description  only  the  mean  of  the  solvated  pro¬ 
tons  and  the  solvation  equilibrium  mean  constant  (for  instance 
a  =  23  =  const.)  can  be  considered  under  isothermal  conditions. 
This  simplification  is  discussed  in  Appendix  A. 

In  terms  of  the  dimensionless  variable 

Ca  Cw  «wt  Cw, 

Va  =  ~p,  <Pw  =  —  X  =  =  —  (2) 

Cwt  Cwt  HA  Ca 


the  local  electroneutrality  condition  and  the  local  hydration  of 
the  membrane  can  be  written  as 


Ca  —  C  A 


dX 

d<Pa  =  ~JJ 


(3) 
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a 

c/Ca  -f-  Cw  =  Cwt  ottpd  -f-  tpw  =  1  d(pw  =  ~  (4) 

kA  dA 

where  a,  <pa,  e,  (pK  can  be  substantially  assimilated  to  the  volu¬ 
metric  fractions  of  the  hydrated  protons  and  water,  respectively. 

The  fluxes  passing  through  the  membrane  depend  on  the  fol¬ 
lowing  driving  forces: 


•  concentration  gradients  inside  the  aqueous  phase,  which 
induce  diffusive  component  fluxes; 

•  voltage  gradients,  which  induce  the  migration  of  the  posi¬ 
tive  ions,  but  are  also  unbalanced  global  forces  acting  on  the 
aqueous  phase; 

•  pressure  gradients,  which  act  on  the  entire  fluid; 

•  hydration  gradients,  which  are  connected  to  capillary  forces 
and  are  then  identifiable  with  an  additional  pressure  gradient 
acting  on  the  entire  fluid. 


By  introducing  an  appropriate  "diffusivity”  for  each  driving 
force,  the  fluxes  of  the  components  (water  and  solvated  protons) 
can  be  expressed  as 


Na  =  ~Dt 


d  Cg 
dz 


„  „  d (VF/RT) 

L'oni'-'oe  , 

dz 


/  DcmCg  \  d  P 
\CmRTj  d^ 


—  (DaxCa) 


dT. 

dz 


(5) 


Nw  =  —D„ 


dCw 

dz 


DvmCv/\  d P  rrt  n  xdl 
CwtRT  )  dz  wk  w  ck 


(6) 


so  that  the  total  water  flux 


Nwi  —  uNa  T  N\\ 


n  ^  d (VF/RT)  DJf  ,nr  dk 

“  nu  dz  RTdz+  lwtdz’ 


(7) 


which  is  a  combination  of  the  preceding  ones,  is  also  consistent 
with  the  momentum  balance  in  the  z  direction.  In  fact,  the  aque¬ 
ous  phase  undergoes  pressure  forces,  capillary  pressure  forces 
and  electrical  forces  acting  on  the  protons,  while  the  electrical 
forces  on  the  anions  are  directly  transferred  to  the  solid  support. 
The  diffusive  forces  balance  themselves  so  that  the  sum  of  the 
diffusive  fluxes  is  zero.  The  other  forces  are  balanced  by  shear 
stresses  under  a  laminar  (or,  better,  a  viscous)  regime,  so  that 
they  can  be  assumed  to  be  proportional  to  the  velocity  and  the 
total  water  flux  Nwt,  as  in  Eq.  (7). 

Moreover,  as  Eq.  (7)  is  a  linear  combination  of  Eqs.  (5)  and 
(6),  the  following  conditions  hold: 


Dac  —  Dwc  —  Dc  I  )- 

Dj]  —  OtDari  D\ 


Ot-DanCa  f  1  a 

Cwt 

otDa\Ca  D-wxC  w 

Cwt 


(8) 


By  introducing  the  dimensionless  variables  and  parameters 
defined  in  the  equations  in  (2)  and  in  the  following: 


NiS 

Cwt  Av 

VF 

~RT 


A.  -  Di 
Cl[  — 

Dc 


it  ■■ 


-  §=  - 

CwtRT  s 


Eqs.  (5)-(7)  can  be  written 

(  1  \  dl  di)  d?r 

na  —  [  da\(pa  T7)  I  ZTZ  (darifa)  , ,  (dan^Pa)^ Z 

v  k-J  d£  df  d£ 


dax  1  \  d  k 

—  ~  y  df 


n  \  d?7 


d 

ak  )  d§ 


dan  \  djT 

r )  d§ 


(9) 


—  —  (10) 


/  a  \  dT.  d7r 

ydy/xtpw  +  ^ J  —  (dW7T(pw)-j^ 


=  - Id- 


Udgi, 

k 


«wt  =  ana  +  «w  =  -di 


d? 

a  \  dl 


dk 


Oldajj;  \ 

7  ) 

d  n 

(11) 

d /; 

drc 

dnd% 

(12) 

wdoijicpa  -f-  dW7T(pw  —  dj{  wdaxcpcL  dvj\(pvj  —  d\ 


d  -  ^ 

uari  — 

a 


(13) 


Some  further  comments  are  useful.  The  description  (5)-(8) 
(or  (10)— (13)  in  dimensionless  form)  of  the  mass  and  charge 
fluxes  passing  through  a  PEMFC  membrane  is  somewhat  more 
detailed  and  more  consistent  than  the  one  presented  and  used  in 
our  previous  works  [6-8].  The  present  description  must  be  con¬ 
sidered,  first  of  all,  as  a  research  tool,  while  the  previous  is  more 
suitable,  given  the  actual  state  of  our  knowledge,  for  applica¬ 
tive  uses,  for  instance  in  simulation.  In  Appendix  B,  however,  it 
will  be  shown  that  Eqs.  (5)-(8)  are  coherent  with  the  previous 
description  and  can  be  reduced  to  it. 

In  particular,  the  fluxes  mentioned  here  refer  to  the  unit  area 
of  the  aqueous  phase,  so  that  some  further  mathematical  steps  are 
necessary  in  order  to  obtain  fluxes  per  unit  area  of  the  membrane 
(see  Appendix  B).  An  analysis  of  the  problem  in  terms  of  the 
transport  properties  of  the  aqueous  phase  has  been  preferred, 
because  this  can  be  considered,  to  some  extent,  independent  of 
the  degree  of  hydration  of  the  membrane.  On  the  other  hand, 
the  hydration  k  must  be  explicitly  considered  in  the  capillary 
forces  and  directly  affects  the  molar  fractions  (pa  (or  volumetric 
fractions  atpa)  of  the  solvated  protons. 

The  consideration  of  a  unique  proton  species  with  a  constant, 
mean,  solvation  number  a  would  seem  to  be  an  over  simpli¬ 
fication.  As  a  matter  of  fact  (see  Appendix  A),  this  choice  is 
completely  correct  only  when  a  thermodynamic  equilibrium 
condition  acts  on  the  proton  species  locally. 

The  generalised  Nernst-Planck  expression,  according  to 
which  the  phenomenological  coefficients  of  the  electrical,  diffu¬ 
sive  and  pressure  forces  have  been  determined,  could  also  be  too 
simple.  For  this  reason,  it  is  worth  noting  that  other  conditions 
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acting  on  the  various  mobility  coefficients,  such  as  those  below, 
could  be  introduced: 

dr,  =  d„=  1  dx^\  (14) 

Roughly  speaking,  the  electrical  forces  per  unit  volume  are 
Cad(VF)!dz  and  the  corresponding  pressure  forces  are  dP/dz, 
while  the  capillarity  forces  are  only  proportional  to  RT dX/dz. 
Therefore,  Dn  and  Dn  could  be  assumed  to  be  equal,  while  D, 
is  certainly  independent.  As  the  diffusive  forces  of  an  ideal  sys¬ 
tem  are  RT dCw/dz,  Dc  could  also  be  assimilated,  at  least  as  a 
first  approximation,  to  Dtl  and  Dn,  as  is  suggested  by  the  ideal 
form  of  the  Nernst-Planck  law.  Nevertheless,  to  avoid  oversim¬ 
plifications,  the  following  discussion  will  be  developed  without 
taking  into  account  the  first  of  the  equations  in  ( 14),  which  will 
only  occasionally  be  referred  to  and  discussed  in  terms  of  orders 
of  magnitude. 

3.  Open-circuit  membrane 


In  the  system  (10)— (13)  the  hydration,  pressure  and  voltage 
gradients  are  usually  considered  as  independent  variables  on 
which  the  fluxes  depend.  On  the  contrary,  the  fluxes  can  also 
be  treated  as  independent  variables  in  order  to  determine  what 
gradients  are  consistent  with  a  given  set  of  fluxes.  In  particular, 
beside  the  equilibrium  condition  of  the  membrane,  in  which  all 
fluxes  as  well  as  all  gradients  vanish  and  the  membrane  voltage 
ought  to  obey  the  Nernst  law,  stationary  conditions,  in  which 
only  one  or  some  fluxes  are  zero,  are  also  possible. 

A  simple  but  particularly  interesting  condition  is  that  of  the 
open-circuit  ( na  =  0),  in  which  Eqs.  ( 10)— (13)  become 

dr]  [(ocdax  -  a/X)(dX/d%)  +  (ada7l)dTX /  d%] 

—  = - ,  na  =  \) 

d£  dr, 

(15) 


/lwt  — 


na  —  0 


adax  o:\ 
1  +  X2) 


dX 

df 


d7T 

df’ 

(16) 


The  open-circuit  condition  does  not  necessarily  correspond  to 
null  gradients,  but  only  imposes  a  relationship  on  the  gradients 
and  the  total  water  flux.  In  other  words,  if  the  membrane  is 
subjected  to  external  driving  forces  in  terms  of  hydration  and/or 
pressure  gradients,  some  related  voltage  gradients  could  also 
subsist  under  open-circuit  conditions.  Pressure  forces,  as  well 
as  capillary  forces  induce  a  water  flux  and,  then,  tend  to  induce 
a  proton-drag  flux;  as  the  latter  is  not  possible  because  of  the 
open-circuit  constraint,  an  antagonist  voltage  gradient  has  to 
occur  as  well. 

The  condition  requiring  that  all  gradients  are  zero  (hydration, 
voltage,  pressure,  temperature),  and  not  that  requiring  only  that 
the  current  is  zero,  is  to  be  considered  sufficient  for  thermody¬ 
namic  equilibrium  throughout  the  membrane.  The  null-current 
(open-circuit)  condition  is  only  necessary  for  the  equilibrium 
and  it  is  consistent  with  a  dissipative  steady-state  (or  a  dissipa¬ 
tive  structure),  according  to  Prigogine’s  terminology  [11]. 


In  brief,  the  immobilisation  of  the  charges  under  the  open- 
circuit  condition  can  co-exist  with  a  water  flux  through  the 
membrane  and  a  corresponding  over  potential.  For  instance, 
under  open-circuit  and  isobar  conditions,  Eq.  (16) 


dr]  _  [(adax  -  a/X)/(dx  -  adak/X  +  a/X2)\ 

df  d,i 


drj 


na  =  0; 
(17) 


indicates  that  a  voltage  loss  through  the  membrane  may  be  asso¬ 
ciated  with  a  negative  water  flux  (that  is  from  the  cathode  to  the 
anode). 

As  the  membrane  is  to  some  extent  also  permeable  to  molec¬ 
ular  hydrogen,  a  slight  positive  hydrogen  flux  (from  the  anode  to 
the  cathode)  can  be  assumed;  this  flux  should  undergo  oxidation 
at  the  cathode  with  the  formation  of  water  and  should  so  induce 
a  counter-flux  of  water  and  an  open-circuit  voltage  loss. 

At  this  point  the  question  is  whether  this  phenomenon  is  truly 
important  or  not;  here  the  discussion  can  only  be  initiated.  In  the 
first  place,  the  quantity  of  open-circuit  water  flux  that  could  cause 
the  observed  voltage  losses  depends  on  the  values  of  the  mobility 
coefficients.  By  allowing  for  the  second  equation  in  (13)  and  the 
assumptions  X  )$>  a  >  1 ,  dax,  dax  of  the  order  of  unity  or  greater, 
Eq.  (17)  can  be  written 


fd ij\  dn  [(X2/a)dwx  ~  Mvi  +  1] _  /  d^N  dn  dwk 

V  d§  )  0  A  Xdax  -  1  V  d£  /  0  «  daX 

(18) 


Then,  a  comparison  with  a  reference  condition  related  to  the 
usual  running  of  the  cell  under  uniform  pressure  and  hydration 


yields 

(^wt)o  ^  7.  (dr]/dt;)0  dwX  ^  ^ 

(«a)r  OC  (dr)/d^\  dal  \daXJ 


(19) 


(20) 


By  assuming  that  (d;//df)r  is  about  three  times  greater  than 
(dri/d%)0  (AVo»0.2V;  AEr«0.6V)  and  X/a^  9,  Eq.  (20) 
demonstrates  that  the  open-circuit  water  flux  responsible  for  the 
observed  voltage  losses  is  very  low  (|(nwt)0/(/7Q.)r|  <<£  1)  only  for 
dwx/dax  1.  On  the  contrary,  if  dwx/dax  ~  1,  Eq.  (14)  demon¬ 
strates  that  the  open-circuit  water  flux  and  the  corresponding 
hydrogen  consumption  should  be  high  enough  to  be  evident. 

Similar  considerations  can  be  made  directly  from  Eq.  (17). 
The  mobility  coefficient  of  the  protons  in  relation  to  the  hydra¬ 
tion  gradients  (du/  )  could  differ  greatly  from  dwX,  so  that  the 
denominator  of  the  second  term  in  (17)  could  be  very  low  and 
the  voltage  loss  could  become  sensitive  to  the  water  flux.  In  this 
case  open-circuit  voltage  losses  of  the  order  of  the  observed  ones 
could  be  foreseen  in  association  with  molecular  hydrogen  and 
water  fluxes  that  are  substantially  low  enough  to  pass  unnoticed. 

The  existence  of  transport  properties  of  the  charged  species 
inside  the  membrane  system  very  different  to  those  of  the  neutral 
ones  is  not  unlikely.  A  hydration  gradient  is  associated  with  the 
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spatial  asymmetry  of  the  porous  system  through  the  membrane 
thickness  (for  instance  the  pore  diameter  will  increase  with  the 
degree  of  hydration),  so  that  asymmetry  in  the  geometric  dis¬ 
position  of  the  negative  charges  fixed  on  the  pore  surface  could 
also  occur.  The  electrical  forces  arising  from  this  asymmetrical 
arrangement  are  likely  to  be  much  greater  than  the  capillarity 
forces,  so  that  the  additional  electrical  forces  due  to  hydration 
gradients  acting  on  protons  (that  is  the  principal  contribution  to 
da/)  could  be  much  greater  than  the  capillary  pressure  forces 
acting  on  neutral  water  (that  is  dw/). 

Secondly,  at  first  sight,  the  direct  consideration  of  the  extent  of 
the  hydrogen  fluxes  also  indicates  that  this  flux  is  rather  unlikely 
to  be  the  only  one  responsible  for  the  relevant  OCVL,  unless  Eq. 
(14)  is  not  true  at  all.  An  estimation  of  the  order  of  magnitude 
of  the  flux  of  the  molecular  hydrogen  through  the  membrane 
and  the  corresponding  counter-flux  of  water  is  presented  in  the 
following  section.  The  hydrogen  flux  is  very  low,  depending  on 
its  solubility  in  the  membrane  system.  If  this  is  assimilated  to  its 
solubility  in  water,  Eq.  (14)  is  assumed  to  hold  and  the  reference 
values  a  =  23  and  X  =  20  are  used.  Eq.  (17)  foresees  voltage 
losses  much  lower  than  0.2  V.  However,  besides  the  question  of 
the  correctness  of  Eq.  (14),  hydrogen  solubility  could  also  be 
strongly  affected  by  adsorption  phenomena  or,  at  least,  by  the 
strong  non-ideality  of  the  aqueous  phase  inside  the  membrane, 
so  that  much  greater  hydrogen  fluxes  could  be  possible. 

In  our  opinion  the  question  of  hydrogen  solubility,  as  well  as 
that  related  to  the  asymmetry  effects  of  the  hydration,  are  worth 
further  theoretical  and  experimental  attention,  as  they  could  lead 
to  a  better  determination  of  at  least  a  part  of  the  open-circuit 
irreversibility. 

4.  Membrane  stability 

When  discussing  fluxes  which  cross  the  membrane  under 
open-circuit  conditions,  another  set  of  questions  appears  to  be 
of  some  interest  and  not  completely  marginal  or  fanciful.  These 
questions  involve  membrane  stability  and  the  origin  of  the  driv¬ 
ing  forces  that  act  on  the  membrane.  One  could  ask  oneself,  in 
particular,  whether  the  absence  of  gradients  imposed  by  external 
constraints  is  equivalent  to  that  of  null  gradients  and  is,  there¬ 
fore,  a  sufficient  condition  for  equilibrium  or,  on  the  contrary,  the 
spontaneous  creation  of  gradients  inside  the  membrane  is  pos¬ 
sible  as  a  result  of  local  instability,  perhaps  in  connection  with 
apparently  secondary  phenomena.  For  instance,  when  recalling 
that  the  charge  vectors  are  at  least  two  (bi-  and  trisolvated  pro¬ 
tons),  the  open-circuit  condition  appears  to  be  compatible  with 
fluxes  in  the  opposite  direction  to  the  two  vectors,  so  that  the 
problem  of  the  stability  of  the  equilibrium  becomes  more  sub¬ 
jective  and  delicate. 

Quantitative  checks  of  thermal  fluxes  (both  conductive  and 
enthalpic),  indicate  that  heat  transfer  through  the  membrane  is 
very  likely  to  be  stable,  which  means  that  the  total  equivalent 
conductivity  of  the  membrane  system  is  positive.  Similar  results 
can  be  obtained,  at  least  at  the  first  approach,  when  the  water 
transfer  is  considered,  that  is  the  total  equivalent  diffusivity  of 
the  water  is  positive  too.  Nevertheless,  it  remains  difficult  to 
definitely  exclude  the  existence  of  any  sensitivity  or  instabil¬ 


ity  phenomenon  able  to  amplify  the  effects  of  the  unbalancing 
of  other,  non-electrical,  driving  forces  in  terms  of  appreciable 
voltage  differences.  The  concentration  gradients  of  molecular 
hydrogen  (or  oxygen)  and  the  affinity  of  the  oxidation  reaction 
are  examples  of  these  kinds  of  forces. 

An  example  of  a  possible  source  of  instability  can  be  given 
as  follows.  In  Eq.  (20),  the  electrical  forces  associated  with 
hydration  gradients,  besides  having  a  greater  absolute  value  than 
the  capillary  ones,  could  even  have  the  opposite  sign,  that  is 
\dwX/dal\  <SC  1  with^wx  <  0.  In  such  a  case,  the  water  flux  asso¬ 
ciated  with  a  voltage  loss  should  go  from  the  anode  to  the  cathode 
and  could  be  related  to  the  occurrence  of  instability  phenomena. 
In  other  words,  if 

dx  =  adctxVa  +  dwx(pw  ^  ct-f-  <  0  (21) 

the  total  water  flux  should  tend  to  go  towards  the  increasing 
humidity,  due  to  the  drag  effect  of  the  solvated  molecules,  and 
should  so  create  a  spontaneous  increment  in  the  hydration  gra¬ 
dient. 

5.  Fluxes  of  neutral  hydrogen 

In  the  last  section  we  have  possibly  tended  to  speculate  too 
much.  To  return  to  more  concrete  arguments,  it  seems  to  be 
useful  to  evaluate  the  order  of  magnitude  of  the  fluxes  of  the 
neutral  reagents,  hydrogen  and  oxygen,  which  cross  a  PEMFC 
membrane  and  cause 

•  the  consumption  of  reagents  without  a  corresponding  produc¬ 
tion  of  energy; 

•  voltage  losses  that  penalise  the  conversion  efficiency  of  the 
overall  electrochemical  reaction. 

These  voltage  losses,  which  also  exist  in  an  open  circuit 
[9,10],  can  be 

•  in  part  closely  associated  with  the  fluxes  of  the  inert  gas  and 
then  be  localised  inside  the  membrane; 

•  in  part  the  result  of  the  interaction  of  inert  fluxes  with  an 
electrode  reaction,  and  so  localised  only  at  the  electrode. 

In  the  following,  the  solubility  of  molecular  hydrogen  and 
oxygen  in  water  will  be  used;  the  correctness  of  all  the  consider¬ 
ations  made  depends  on  the  accuracy  of  such  data  in  representing 
the  solubility  of  the  membrane  system.  The  solubility  of  hydro¬ 
gen  in  water  at  temperatures  similar  to  those  occurring  in  a 
cell  is  very  low.  Henry’s  constants  for  hydrogen  and  oxygen 
at  100  °C  are,  respectively,  7.45  x  104  and  7.01  x  104atm.  At 
an  atmospheric  PEMFC  anode,  reference  can  be  made  to  about 
1  atm  of  hydrogen,  with  an  equilibrium  molar  fraction  of  about 
1.3  x  10-5,  while  at  the  cathode  the  oxygen  partial  pressure  is  of 
the  order  of  0.2  atm  or  lower,  with  an  equilibrium  molar  fraction 
lower  than  2.8  x  10-6.  Even  allowing  for  the  stoichiometry  of 
the  oxidation  reaction,  hydrogen  exceeds  oxygen  and  it  is  there¬ 
fore  reasonable  to  make  reference  to  a  molecular  hydrogen  flux 
which  dissolves  at  the  anode,  migrates  through  the  water  phase 
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inside  the  membrane  and  ends  up  reacting  with  the  oxygen  at  (or 
near)  the  cathode,  so  producing  water.  This  water  can  counter- 
diffuse  towards  the  anode  and  cause  a  (little  ?)  voltage  loss  in 
the  membrane. 

Due  to  its  very  low  solubility,  this  molecular  hydrogen  flux 
must,  in  turn,  be  very  low.  Making  reference  to  an  effec¬ 
tive  diffusivity  of  the  membrane  of  10-6  cm2  s-1  and  a  mem¬ 
brane  thickness  of  0. 1  cm  a  hydrogen  flux  N'Hi  of  the  order  of 
7  x  10-12  mol  cm-2  s_1  is  obtained,  which  corresponds  to  a  cur¬ 
rent  of  the  order  of  1.4  x  10-6  Acm~2. 

The  insignificance  of  these  neutral  hydrogen  fluxes  in  a  closed 
circuit  should  make  it  possible  to  ignore  their  cross-effects  on 
the  other,  much  more  important  fluxes  of  water  and  protons.  On 
the  other  hand,  the  persistence  of  such  irreversible  phenomena 
in  an  open  circuit  makes  it  necessary  to  investigate  the  matter 
further. 

As  already  said  in  the  previous  section,  hydrogen  fluxes  of 
such  an  order  of  magnitude  could  induce  significant  voltage 
losses  through  the  membrane  only  in  association  with  impor¬ 
tant  asymmetries  in  the  conductivity  coefficients  and  the  pore 
geometry.  It  still  remains  to  consider  what  happens  outside  the 
membrane  and,  in  particular,  what  effects  could  be  induced  on 
the  cathode,  by  a  hydrogen  flux  and  its  consequent  oxidation 
reaction. 


current  Ioc  =  1.4  x  10-6  A  cm-2,  an  exchange  current  value 
70  =  3.5  x  10-9  A  cm-2  should  be  obtained.  Therefore,  a  typi¬ 
cal  PEMFC  cathode  could  be  described  in  one  of  the  following 
ways. 

(a)  The  cathode  is  characterised  by  an  exchange  current 
7*  =  /oc  =  1.4  x  10-6  Acm-2 
and  an  open-circuit  voltage  loss 
Vo  =  0.2  V 

which  is  attributed  to  imprecisely-stated  dissipative  phe¬ 
nomena.  Further  voltage  losses  under  closed-circuit  con¬ 
ditions  follow  Tafel’s  law: 


(b)  The  cathode  is  characterised  by  a  much  lower  exchange 
current 

70  =  3.5  x  10-9  Acm-2 


6.  Open-circuit  electrode 

The  fact  that  such  low  fluxes  cause  such  significant  irre¬ 
versibility  effects  suggests  that  further  study  should  also,  or 
prevailingly,  be  oriented  towards  strongly  non-linear  and  sen¬ 
sitive  phenomena,  such  as  the  activation  mechanisms  on  the 
electrodes. 

The  neutral  oxygen  and  hydrogen  fluxes  that  enter  the  cathode 
and  undergo  an  oxidation  reaction  there  correspond  to  a  parasitic 
current  completely  inside  the  cathode.  Then,  in  a  closed  circuit, 
the  net  current  can  be  considered  as  the  difference  between  the 
total  current  7t  and  the  parasitic  current  70C 

I  =  h~Io  c  =  Mex  p[fi(v  +  Vo)]  ~  exp[-(yS  -  1)0;  +  t?0)]} 


and  a  molecular  hydrogen  flux  corresponding  to  a  very  low 
parasitic  current 

7oc  =  1.4  x  10-6  Acm-2 

with  which  a  voltage  loss  of 

Vo  =  0.2  V 


is  associated.  The  total  open-circuit  voltage  losses  follow 
Tafel’s  law: 


V  +  Vo  =  In 


with  v  =  In 


and 


-70{exp[^0]  -  exp[-(/3  -  l)ij0]}  (22) 

The  parasitic  current  is,  in  turn,  connected  to  a  voltage  loss 
Vo ,  which  also  exists  in  an  open  circuit,  when  the  net  current 
vanishes,  together  with  the  additional  voltage  loss  v' 

Ioc  =  70{ex p[/3/70]  -  exp[-(yS  -  l)r;0]}  (23) 

By  assuming  /l  =  0.5  and  Vo  ~  12,  that  is  an  OCV  voltage 
loss  of  about  0.2  V,  the  preceding  equations  can  be  reduced  to 
the  Tafel  form 

Ioc  =  Io  exp(/S?;0)  (24) 

I  =  70  exp(fivo)[e*p(Pv)  ~  1]  =  /0exp(/3/;0)exp(/3(;) 

=  Ioc  exp(/b;)  (25) 

By  making  reference  to  the  aforesaid  indicative  value  of 
the  molecular  hydrogen  flux  and  the  corresponding  parasitic 


The  two  descriptions  are  identical  and  differ  only  in  the  terms 
used,  which  seem  to  be  more  analytical  and  interpretative  in  the 
second,  in  which  the  possibility  of  identifying  the  irreversible 
electrode  phenomena  under  open-circuit  conditions  in  terms  of 
a  “parasitic  current”  is  proposed. 

The  important  point  is  that  this  parasitic  current  is  here  asso¬ 
ciated  with  a  hydrogen  flux  that  crosses  the  membrane,  under 
steady-state,  non-equilibrium,  open-circuit  conditions.  On  the 
other  hand,  only  more  accurate  experimental  investigations, 
aimed  at  measuring  small  open-circuit  hydrogen  fluxes  and  inter¬ 
nal  currents  would  make  it  possible  to  verify  or  refute  this 
interpretation.  While  we  are  waiting  for  more  concrete  proof, 
the  use  of  the  term  “parasitic  current”  to  explain  the  open-circuit 
voltage  losses,  or  a  part  of  them,  remains  as  vague  as  the  appar¬ 
ently  more  generic  “dissipative  phenomena”. 
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7.  Conclusions 

The  performance  of  polymeric  electrolyte  membrane  fuel 
cells  strongly  depends  on  transport  mechanisms  and  dissipa¬ 
tive  phenomena  that  occur  inside  the  membrane  and  the  elec¬ 
trodes.  This  work  could  represent  an  interpretative  contribution 
to  understanding  open-circuit  voltage  losses,  which  can  be  up 
to  30%  of  the  effective  voltage  of  an  operating  cell  and  cause 
considerable  damage  from  a  technological  point  of  view. 

In  particular,  the  transport  mechanisms  taking  place  inside 
the  membrane  and  the  interaction  of  these  with  the  adjacent 
electrodes  have  been  analysed  in  detail  to  better  understand  and 
control  cell  performance. 

In  an  open  circuit,  low  hydrogen  fluxes  could  be  associated 
with  strongly  non-linear  and  sensitive  phenomena,  such  as  the 
activation  mechanisms  on  the  electrodes.  An  analytical  and  inter¬ 
pretative  description  has  been  proposed  to  try  to  identify  the 
irreversible  electrode  phenomena  under  open-circuit  conditions 
in  terms  of  parasitic  currents. 

Therefore,  we  assume  that  parasitic  currents  are  linked  to  a 
hydrogen  flux  which  would  cross  the  membrane,  under  steady- 
state,  non-equilibrium,  open-circuit  conditions:  only  very  accu¬ 
rate  experimental  investigations  aimed  at  measuring  small  open- 
circuit  hydrogen  fluxes  and  internal  currents  would  make  it 
possible  to  verify  the  interpretation  proposed  in  this  work. 
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and  the  definition  of  the  mean  degree  of  hydration  (or  solvation) 
of  the  protons  can  be  stated  as 
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considered  here  as  a  function  of  the  local  hydration  of  the  mem¬ 
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or,  in  dimensionless  terms, 
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Appendix  A 


If  the  bihydrated  and  trihydrated  protons  are  considered  sepa¬ 
rately  and  the  ratio  of  their  respective  concentrations  is  assumed 
to  obey  an  equilibrium  relationship  of  the  type 

C3  =  k(X)C2  (A.l) 


instead  of  that  in  Eqs.  (3)  and  (4),  the  following  are  obtained: 


1 

<P2  +  (03  =  (02(1  +  k)  =  - 

A 

dk 

d(p3  =  -(1  -  b)- 2 


dq)2  =  -b 


dk 

k2 


(A. 2) 


2(p2  +  3(03  +  <Pw  —  <P2(2  +  3k)  +  (0W  =  1 
dk 

d(0w  =  (3  -  b)- 2 


where 


b  = 


1  +  k  +  X  dk/dX 
(1  +k)2 


(A. 3) 


(A.4) 


b  \  dk  dr]  d7T 
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d? 

(A.  10) 

3-b 

dk 

dn 

nw  — 

H"  “wc  ^2 

—  (dwjziPvj)—— 

d§ 

(A.  11) 

Obviously,  by  adding  Eqs.  (A.9)  and  (A.  10),  Eq.  (10)  is  once 
again 

'da X  _  J_\  dk  _  / donj\  dj7  _  / dom\  dn 

k  X7-)  d§  U/d?  U  J  d$ 

under  the  conditions 

bd2C  +  (1  -  b)d3c  =  1  d2x  +  kd3x  =  (1  +  k)dax 


dx,]  T  kd3 tj  —  (1  T  k)dc/yj  d2n  T  kd3^  —  (1  T  k )rla7T 


(A.  12) 


Similarly,  by  combining  Eqs.  (A.9)-(A.l  1)  in  the  total  water 
flux 


nwt  =  2n2  +  3«3  +  nw 
Eq.  (12)  is  once  again 

dk  1  d??  d7r 

n wt  =  an a  +  nw  =  -dx  ~rr  —  tt  ~tt 
d^  k  d^  d§ 


(A.  13) 


(12') 
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under  the  further  conditions 
2 bd2c  +  3(1  -  b)d3c  -  (3  -  b)  =  0 


2 d2r]  +  3kd2tl  =  (1  +  k)  d^ 

2d2n  +  3kd2n  —  (2  +  3k)dVJ7I  =  X(d7r  —  rfwjr)(l  +  k ) 


2d2x  +  3kd2x  -  (2  +  3k)dwx  =  X(dx  -  dwx)(l  +  k) 

The  conditions  (A.  12)  and  (A.  14)  show  that,  as  long  as  the 
local  equilibrium  condition  (A.  1)  holds,  univocal  links  also  hold 
among  the  mobilities  of  the  proton  species  (for  instance  d2>1  and 
d2v)  and  that  of  the  mean  solvated  species  (for  instance  dri). 

Appendix  B 


Eqs.  (5)  and  (6)  aim  to  describe  mass  and  charge  transport 
phenomena  occurring  in  the  aqueous  phase  of  a  heterogeneous 
membrane  system.  The  fluxes  and  the  concentrations  refer  to  the 
aqueous  phase,  so  that  A;  is  moles  of  i  per  unit  of  aqueous-phase 
area  and  C,  is  moles  of  i  per  unit  of  aqueous-phase  volume. 
Moreover,  in  Eqs.  (5)  and  (6)  the  gradients  dC,/dz  (that  is  the 
diffusive  forces)  and  dX/dz  (that  is  the  capillary  forces)  are  con¬ 
sidered  separately.  So,  two  steps  have  to  be  made  in  order  to  bring 
the  present  description  back  to  the  one  used  more  pragmatically 
in  our  simulation  works. 

As  a  first  step,  by  using  Eqs.  (2)-(4),  the  concentration  gra¬ 
dients  and  the  hydration  gradients  can  be  transformed  into  each 
other  by  the  differential  relationships 

dl  =  -  (1)  dc“  =  (st) dc"  (B1) 


So,  by  using  the  equations  in  (B.l),  Eqs.  (5)  and  (6)  can  be 
written 


Na 


dCc, 

dz 


F  \ 

Dar,Ca—  J 


dV 

dz 


Ca  \ 

c^rt) 


d  P 
dz 


(B.2) 


Aw 


DwxCwa<%)  J 


dCw 

dz 


/  AyttCw  \  dE 

\CwtRTj  dz 

(B.3) 


In  the  second  step  the  fluxes  A,  and  the  concentrations  C; 
have  to  be  transformed  into  the  fluxes  N[  and  the  concentrations 
C'r  which  refer  to  the  whole  heterogeneous  membrane  system 

NiS  =  a;  CiS  =  C\  (B.4) 

by  using  the  volumetric  fraction  s  of  the  aqueous  phase  in  rela¬ 
tionship  to  the  volume  of  the  membrane  system.  This  volumetric 


fraction  is  directly  linked  to  (and  roughly  coincides  with)  the 
swelling  parameter  e 


s  = 


1  +  e 


(B.5) 


and  can  also  be  expressed  in  terms  of  the  concentration  C'a 

1  -  CL  1  -CL 


s  = 


C' 


C' 


(B.6) 


As 


so  that  the  following  differential  relationships  between  dC,  and 
d C':  can  be  obtained: 


d CL  =  s  dCa 


d  CL  = 


s  4~  hs- 


1  —  e 


dCw 


(B.7) 


By  substituting  Eqs.  (B.4)  and  (B.7)  in  (B.2)  and  (B.3)  and 
recalling  (2)-(4),  the  final  form 


A  = 


Dgc  kDgX 


Dgjre\  d P 


dC 

dz 


F  \dV 

Og,Cwts— j  — 


XRT  J  dz 


(B.8) 


A  = 


OlD-wc  ~\~  OL)h 

dc; 

d  P 

X  —  (k  —  a)s 

dz 

XRT 

dz 

(B.9) 


of  the  flux  equations  is  obtained,  which  has  the  same  structure 
as  Eqs.  (10)  and  (11)  of  the  reference  [6]  and  can  easily  be 
compared  to  them.  The  transformation  of  one  set  of  diffusivity 
coefficients  into  the  other  is  now  quite  straightforward. 
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